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We present an extensive numerical study of the ferromagnetic Kondo lattice model with quantum 
mechanical 5 = 3/2 core spins. We treat one orbital per site in one dimension using the density 
matrix renormalization group and include on-site Coulomb repulsion between the electrons. We 
examine parameters relevant to manganites, treating the range of low to intermediate doping, < 
X < 0.5. In particular, we investigate whether quantum fluctuations favor phase separation over the 
ferromagnetic polarons observed in a model with classical core spins. We obtain very good agreement 
of the quantum model with previous results for the classical model, finding separated polarons which 
are repulsive at short distance for finite t2g superexchange J'. Taking on-site Coulomb repulsion 
into account, we observe phase separation for small but finite superexchange J', while for larger J' 
polarons are favored in accordance with simple energy considerations previously applied to classical 
spins. We discuss the interpretation of compressibilities and present a phase diagram with respect 
to doping and the t2g superexchange parameter J' with and without Coulomb repulsion. 
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I. INTRODUCTION 



The ferromagnetic Kondo lattice model has been 
widely used as a minimal model to describe some features 
of the manganitesi Lai_2:Sr2:Mn03, Lai_a;Sri+a;Mn04, 
and La2-2a;Sri_|_2a;Mn207. The model contains one itin- 
erant Cg orbital and the t2g core spin at each site. The eg 
electrons are ferromagnetically coupled to the S — 3/2 
core spins generated by the fully occupied t2g orbitals. 
The large ferromagnetic Hund's rule coupling leads to 
the formation of two bands, the lower and the upper 
Kondo band, with the Cg electrons predominantly par- 
allel to the core spins in the lower band. The core spins 
strongly influence the mobility of the Cg electrons via 
double exchange (DE). At high hole density, this leads 
to a ferromagnetic arrangement of the core spins, while 
antiferromagnetism is preferred for the completely filled 
lower Kondo band. In the opposite case of an empty 
Kondo band, the t2g core spins are antiferromagnetically 
oriented due to superexchange. 

Since the core spins have a fairly large spin, S — 
3/2, they are frequently approximated by classical spins, 
greatly simplifying calculations. Furthermore, the on-site 
Coulomb repulsion U between the eg electrons is often ne- 
glected because double occupancy is already suppressed 
by the Hund's rule coupling and because its treatment 
considerably increases the numerical effort;^ A review of 
these semi-classical simulations can be found in Refs.00 
and references therein. Several of these studies^ have 
found that phase separation into regions with either 
ferromagnetically or antiferromagnetically aligned core 



spins occurs when the lower Kondo band is nearly empty 
(n > 0) or nearly filled (n < 1). 

In previous work by some of the current authors, also 
treating the core spins classically, similar numerical re- 
sults were obtained. ^'^ However, a closer analysis of the 
data revealed that the features which had been inter- 
preted by other authors to indicate phase separation (a 
discontinuity in the electron density versus the chemi- 
cal potential, a pseudogap in the one-particle density of 
states) were, in fact, due to small, independent ferromag- 
netic polarons. Likewise, ferromagnetic polarons have 
been found for the almost empty lower Kondo band, for 
localized S = 1/2 quantum spinsi^ For the antiferro- 
magnetic Kondo modelfi^i small ferromagnetic droplets 
were predicted by Kagan et alp^i The question arises 
whether the correct quantum mechanical treatment of 
the 5 = 3/2 core spins would favor phase separation in- 
stead of independent ferromagnetic polarons, especially 
for r = 0. In this paper, we address the influence of 
quantum spins on this issue. 

The impact of a quantum mechanical treatment of 
spins on models for the Manganites has been addressed 
in Refs . IT^ITMT^ITsl . In one dimension, quantum mechan- 
ical core spins with 5 = 1/2 were employed in a number 
of studies conducted with the density matrix renormal- 
ization group (DMRG)i^*2iiSiii Recently, Garcia et ali^ 
presented a phase diagram for S = 1/2, which, however, 
was determined for only three values of the density and 
did not address the physically interesting region of dop- 
ing, X < 1/3, treated in this paper. Quantum mechanical 
spins with 5 = 3/2 have been briefly addressed in an ex- 
ploratory studyi^ The authors report phase separation 
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when the lower Kondo band is nearly filled. 

In this work, we present extensive calculations for the 
one-dimensional (ID) ferromagnetic Kondo lattice model 
with S — 3/2 core spins using the DMRG. We observe 
that quantum spins yield results in very good agreement 
with previous calculations^ for a model with classical core 
spins: for finite t2g superexchange J', polarons are fa- 
vored over phase separation. In addition, we include an 
on-site Coulomb repulsion U between the Cg electrons: 
accordingly, the superexchange parameter J' which fa- 
vors antiferromagnetic alignment of the core spins is 
renormalized and the effective antiferromagnetic coupling 
Jeff is weakened. This leads to an increased polaron size 
for large J' and phase separation for small J'. 

The remainder of this paper is organized as follows. 
In Sec. ^ we define the model Hamiltonian, for which 
we present DMRG results in Sec. IIIII We discuss the 
ground state configurations in Sec. IIII Al In Sec. IIII A II 
we show that the polarons are actually repulsive at short 
distance. We then discuss the impact of the antiferro- 
magnetic t2g superexchange J' and of the Hubbard repul- 
sion U and show that quantum mechanical core spins are 
very well-approximated by classical spins fSec. IIII A 2|l . 
We discuss the transition to the homogeneous ferromag- 
netic chain in Sec. IIII A3l and present a phase diagram in 
Sec. IIII B1 Negative compressibility and the discontinuity 
in the density vs. the chemical potential arc often taken 
as a sign for phase separation. In Sec. IIII CI we argue 
that negative compressibility is an uncertain result when 
obtained from numerical methods such as the DMRG, 
quantum Monte Carlo, or exact diagonalization and show 
that a discontinuity in the density can equally well result 
from small independent polarons. Finally, Sec. IIVI sum- 
marizes and discusses the results presented in this work. 



Hund's rule coupling between the Cg electrons and the 
t2g core spins. In this work, we take Jh = 8, which 
corresponds to Jh = 6 in Refs. and 0, if one com- 
pensates for the normalization of classical core spins to 
iS*! = 1. The last term describes an additional direct 
superexchange between the core spins. For manganites, 
this effective interaction favors antiferromagnetic order- 
ing of the core spins, i.e., J' > 0. We vary J' from 
J' = to J' = 0.02. Note that J' = 0.01 corresponds to 
J' = I • 0.01 « 0.02 in the units of Refs. ITgl and^. 

We employ the density matrix renormalization group, 
keeping up to 1000 states at each DMRG iteration and 
treating chains of length L — 24: with open boundary 
conditions. The discarded weight is at most of order 10~^ 
for the results presented here. Our calculations indicate, 
however, that 48 states, as used in Ref. [6|, would lead to 
insufficient accuracy for the system sizes treated here. 



III. DMRG RESULTS 
A. Ground state configurations 

1. Polarons are repulsive 

Figure ^ shows the on-site electron density and the 
nearest-neighbor core spin correlation (S'i5'i+i) for a 
DMRG ground state for a chain of length L = 24 for 
J' = 0.01 and U ~ with three holes. The holes clearly 
form three individual polarons, each extending over ap- 
proximately three sites with ferromagnetically aligned 
core spins that are embedded in an antiferromagnetic 
background. 



II. MODEL HAMILTONIAN AND METHOD 

We study the ferromagnetic Kondo lattice model with 
localized quantum core spins S — 3/2 and one orbital 
per site, including a small Heisenberg-like superexchange 
between the core spins as well as an on-site Coulomb 
repulsion U between the eg electrons: 

~JHY.S,■S, + J'Y,S^■S,, (1) 

i (ij) 

where c^^ (cj^) creates (destroys) an eg electron with spin 
a at site i, fii^cr — c^^r is the corresponding density op- 
erator. Si is the core spin at site i and Si the electron 
spin. The first term describes the electron hopping; the 
hopping integral t — I will be used in the following as the 
unit of energy. The second term describes the Coulomb 
repulsion for the eg electrons; we will treat U = and 
U = 10. The third term describes the ferromagnetic 
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FIG. 1: Local density {ui) and core spin-core spin correlations 
(SiSi+i) of the DMRG ground state for a chain of length 
L = 24 with J' = 0.01 and C/ = 0. Three well-separated 
polarons can be clearly seen. 

Since the DMRG is a variational method, convergence 
to the true ground state is not guaranteed. Therefore, 
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we have checked the consistency of the ground state by 
trying different system buildup strategies with respect 
to particle injection for many parameter sets. Thereby, 
we have examined whether there is degeneracy of states 
with respect to polaron positions and how far the exis- 
tence of well-separated polarons and their positions are 
determined by the details of the DMRG algorithm. In 
addition, one can perform calculations with different z- 
components 5'*°* of the total spin S"°* : increasing 5'*°* 
reduces the size of the Hilbert space. This, in gen- 
eral, allows for more accurate results, while the smaller 
Hilbert space still contains the ground state as long as 

In an L = 24 system with two holes (J' = 0.01, IJ = 
0), we always find two polarons regardless of how the 
system is built up. All energies obtained are degenerate 
within the estimated numerical accuracy. The polaron 
position, however, depends on where the particles are in- 
jected. In all cases, we find polarons that are separated 
by at least a few lattice sites with AFM order. Even 
if we add both holes simultaneously during the system 
buildup, they separate into two polarons. Therefore, we 
conclude that well-separated polarons are effectively in- 
dependent, while they seem to repel each other at short 
distances. 

We estimate the energy connected to this repulsion by 
introducing small electrostatic potentials i?pot = —0.1 in 
order to trap the holes at sites xi and X2. For all dis- 
tances d = \x2 — xi\ > 3, we obtain FM polarons covering 
three lattice sites with on-site densities symmetric with 
respect to xi and X2- The ground-state energy as a func- 
tion of the distance d — \x2~xi \ is shown in Fig. El which 
corroborates the fact that polarons separated by two or 
more sites are effectively independent. The configura- 
tion with only one intermediate site is AE w 0.01 higher 
in energy, and in order to obtain a state with adjacent 
polarons AE « 0.03 has to be paid. 

We therefore conclude that the holes actually have a 
repulsive interaction at short range and that separated 
polarons are energetically favored over phase separation 
into larger FM and AFM regions near the half-filled 
Kondo band for J' = 0.01 and [7 = 0. 

The fact that we can obtain degenerate states with a 
localized polaron at different positions upon doping with 
one hole suggests that polarons are quasi-particles having 
a large effective mass. For physical reasons, we would 
expect that a single polaron delocalizes, forming a heavy 
Fermi liquid, but it appears that polarons have a band- 
width that is too small to be resolved by the DMRG. 



2. Influence of J' and U and comparison with classical 
spins 

In Ref. the energy for phase separation was com- 
pared to that for independent polarons as a function of 
an effective superexchange J^s of classical core spins for 
the almost filled lower Kondo band {n < 1). At small 




Distance d = Ix^ - x^l 



FIG. 2: Ground-state energy Eq as a function of the distance d 
of two polarons on an L = 24-site chain with J' — 0.01, i7 = 0. 
The holes are centered at the position of electrostatic impu- 
rity potentials Bpot = —0.1 and exhibit a symmetric density 
distribution. The energies for d > 5 are degenerate within 
the estimated numerical accuracy as designated by the error 
bars. 
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FIG. 3: (Color online) Dressed core spin correlation function, 
Eq. ||3J , for classical core spins with effective spinless fermions 
for J' = 0, J7 = 0, inverse temperature P — 100 and an 
L — 24-site chain. The ferromagnetic area grows with doping, 
indicating phase separation. 



values of Jcs, an increase in the optimal polaron size was 
found and bipolarons and phase separation dominated. 
The effective exchange Jes is defined in terms of the t2g 
superexchange J' and the energy ecx for virtual excita- 
tions into the upper Kondo band: 



^off ~ J' + . 



(2) 



Without Coulomb repulsion, £cx is given by the energy 
for the low-spin state For finite Hubbard U, the 

energy for the virtual excitations is higher, however, be- 
cause they lead to doubly occupied sites for n « 1; Ecx 
is then given by the energetically higher '^E and 
states.^^ Taking into account Coulomb repulsion should 
therefore have an effect similar to reducing J'. 

The spinless fermion model with classical core spinsiS 
forms polarons for J' = 0.01 (corresponding to J' « 0.02 
in the units of Ref. 0), but for J' = the polarons tend to 
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FIG. 4: (Color online) Dressed core spin correlation, Eq. 0, 
for quantum mechanical core spins for (7 = and (a) J' = 
0.01 and (b) J' = on an L = 24-site chain. In the polaronic 
regime (a) there is no dependence of the size of the ferromag- 
netic area on doping, in contrast to the phase separated case 
(b). 



attract each other and phase separate. This can be seen 
by examining the dressed core spin correlation function 



FIG. 5: (Color online) Dressed core spin correlation, Eq. 
for quantum mechanical core spins for U = 10 and (a) J' = 
0.01 and (b) J' = on an L = 24-site chain. The pronounced 
dependence on doping indicates phase separation. 



fore increases the tendency to phase separation. In fact, 
phase separation takes place even for J' — 0.01 which is 
reflected in the dressed core spin correlations in Fig. \S\ 
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where = (1 — ni|)(l — Uii) gives the hole density 
relative to the half filled chain, i.e., it is only nonzero if 
the site is unoccupied. We have evaluated this observable 
for the classical model with J' = 0. The result, depicted 
in Fig. shows that the ferromagnetic regions around 
the holes grow with doping, i.e., that the polarons attract 
each other and tend to phase separate. 

The dressed core spin correlation function Sh{r) is 
shown for quantum mechanical core spins and U — 
in Fig. 21 For J' = 0.01 the size of the ferromagnetic 
regions around the holes does not grow with doping in 
accordance with our analysis in Sec. IIII jTTl In the case 
of J' = 0, we observe phase separation also for the quan- 
tum model, which is reflected in the pr onoun ced increase 
of FM correlations with doping in Fig. 4(b) Without an 



on-site Coulomb repulsion U, we find that the spinlcss 
fermion model^ agrees qualitatively very well with the 
present DMRG results, confirming that localized S — 3/2 
spins are well-approximated by classical spins. 

If ones takes a closer look at the numerical values of the 
dressed core spin correlation function for both models, 
one has to keep in mind that states with double occupa- 
tion have been projected out in the classical case, while 
such states are included in the DMRG calculations. In 
the latter case, there is a finite hole density (n^) « 0.014 
at sites i within the AFM background far from polarons, 
and still higher hole densities can be observed in the prox- 
imity of polarons. These contributions of the AFM back- 
ground to the dressed core spin correlation function Q 
partly cancel the FM terms from within the polarons in 
the quantum model; this accounts for the almost vanish- 
ing nearest-neighbor correlation Sh.{l)/Sh{0) ~ for one 



hole shown in Fig. 4(a) 

The inclusion of a finite Hubbard U — 10 decreases 
the effective exchange Jeff as discussed above and there- 



3. Transition to the homogeneous ferromagnetic phase 

We first examine the polaronic case J' = 0.01 and 
U — 0, where the polarons extend over approximately 3 
sites. The ground state of an L — 24: chain for x = 1/4 
(6 holes) shown in Fig. 6(a) corresponds to a periodic 



arrangement of polarons, in agreement with the energy 
considerations of Sec. IIII ATI Similar 'island phases' have 
also been found for S — 1/2 core spins at commensurate 
fillings: They consist of small ferromagnetic islands that 
are aligned antiferromagnetically or have one anti-aligned 
spin between them. Adding one more hole, we find that 
the polaronic configuration is destroyed and a large FM 
region forms with two AFM arranged spins at each end, 
see Fig. |6(b)| This could be a sign of phase separation, 
but much larger chains would have to be treated in order 
to clarify this issue. For dopings of a: = 0.375 (9 holes on 
an L = 24 chain) we observe complete FM polarization. 
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FIG. 6: On-site density {ui) and core spin-core spin corre- 
lations (SiSi+i) of the DMRG ground state for J' = 0.01, 
[/ = and (a) 6 holes {x = 1/4) and (b) 7 holes. 

When the core spin supcrexchangc is increased to 
J' — 0.02 (still neglecting Coulomb repulsion U = 0), po- 
laronic states become stabilized up to dopings ol x — 1/3 
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and we find another islan d pha se of AFM stacked po- 
larons for 8 holes; see Fig. 7(a) The result for 9 holes 
in Fig. 7(b) suggests phase separation between an island 
phase and a FM region; however, chains of length L — 24 
are too small to make definitive statements. As the treat- 
ment of much larger systems is not feasible using cur- 
rently available computational resources and the main 
focus of this paper is the almost-filled lower Kondo band, 
we did not investigate the nature of the phase transition 
further. We finally note that upon doping the L = 24 
chain with 10 or more holes, complete FM polarization 
is obtained (not shown). 
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FIG. 7: On-site density (rii) and core spin-core spin corre- 
lations {SiS^+l) of the DMRG ground state for J' = 0.02, 
[/ = and (a) 8 holes (x = 1/3) and (b) 9 holes. 



For parameters J' — 0.02 and U = 10, the polaronic 
regime extents up to a; = 1/4 (6 holes) , so the transition 
to ferromagnetism is very similar to the one obtained 
for J' = 0.01 and U = 0. We find an island phase as 



already shown in Fig. 6(a) and for 7 holes a configuration 
similar to the one shown in Fig. 6(b) At dopings of 8 
holes {x = 1/3) complete FM polarization is obtained for 
J' = 0.02 and U = 10. 



B. Phase diagrams 
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FIG. 8: (Color online) Phase diagram in the n-J' plane for 
U = 0. The symbols designate the following characteristics: 
P: polarons, ILP: island phase (periodic arrangement of po- 
larons), PS: phase separation, FM: ferromagnetic, FR: central 
ferromagnetic region with antiferromagnetically oriented sites 
at the ends of the chain. The symbols (letters) are determined 
for an L = 24-site chain. 



We summarize the impact of the parameters J' and U 
in phase diagrams in the n-J' plane, where n — 1 — x des- 
ignates the filling and J' refers to the t2g superexchange. 
Figure |S1 shows the phase diagram for an L = 24-site 
chain without Coulomb repulsion, U = 0. One sees a po- 
laronic region (P) near the filled lower Kondo band, in- 
cluding periodic arrangement of polarons (island phase, 
denoted as ILP) for commensurate fillings. For large 
J' > 0.01, this region extends to hole densities x fa 1/3, 
while phase separation is found for vanishing J' = 0. 

Some of our results indicate that there might be phase 
separation for intermediate fillings between polaronic and 
FM phases at J' > 0. However, much longer chains would 
be needed to clarify this issue; see also our discussion in 
Sec. 1111 A 31 Corresponding states with a large central 
FM region and eit her so me AFM sites or polarons at 
both ends (cf. Figs. 6(b) and 7(b) I have been designated 
as "FR" in the phase diagram. Regions are labeled "FM" 
if all nearest-neighbor spin correlations are positive. 

For U — and J' = 0, we find phase separation (PS): 
By varying the system buildup, we are able to obtain 
both polaronic configurations and states with a larger FM 
region embedded in an AFM background; for these pa- 
rameters, the phase separated states always have a lower 
energy. 

While we expect macroscopic polarization in the ther- 
modynamic limit near the empty Kondo band, as also 
observed in ID for the t — t' — U model22i2^ and a Cu-0 
chain^^, we did not investigate the extent of polariza- 
tion for finite systems in detail. We note that saturation 
should not be expected for quantum spins in the thermo- 
dynamic limit according to Ref. ITsl 
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FIG. 9; (Color online) Phase diagram as in Fig. m but with 
[/= 10. 

Fig.inishows the phase diagram for Coulomb repulsion 
U — 10. Comparing with U — Q, we see that the po- 
laronic phase is suppressed at lower doping and phase 
separation takes place even for J' = 0.01. This is con- 
sistent with the reduction in the effective parameter Jeff 
with increasing U . For parameters J' = 0.01 and U = 10, 
we find that two holes can gain Ai? w 0.008 forming a 
bipolaron instead of two separated polarons. Compared 
to the energy E « —137.318 of the bipolaronic state, the 
energy difference is rather small and this is one reason 
for the sensitivity of the DMRG with respect to system 
buildup: If there are states very close in energy but far 
apart in phase space, numerical results will strongly de- 
pend on the initial configuration. We have designated 



6 



the corresponding parts with "PS" in Fig. |U| 



C. Compressibility 

The inverse compressibihty of a system can be com- 
puted approximately by numerical differentiation of en- 
ergies by 

_i _ EjNe + A, L) + EjN, - A, L) - 2E{N,, L) 

" ~ ~ A^ ' 

(4) 

where E{Nf.,L) is the total energy of a chain with N^, 
electrons on L sites and A is the difference in particle 
number. Negative values are sometimes taken to be an 
indication of phase separation>& 

We want to argue that this condition is neither nec- 
essary nor sufficient: On the one hand, negative values 
only result from finite-size effects and should vanish in 
the thermodynamic limit leading to 0, as will 

be discussed below. On the other hand, we will show in 
Fig.^lthat polarons can cause k^^ « as well. Such ob- 
servations therefore have always to be complemented by 
an investigation of correlation fimctions, e.g., the dressed 
core-spin correlation, Eq. (PJ, in order to show the exis- 
tence of two distinct phases. 

First, we note that numerical methods like the DMRG 
or Monte Carlo calculations which do not impose homo- 
geneity, as for instance mean-field theory would, should 
yield separated phases in different spatial regions in a 
proportion that minimizes the free energy — in effect, 
the system performs the Maxwell construction by itself 
and should thus avoid negative compressibilities in the 
thermodynamic limit. 

For phase separation on finite systems, however, neg- 
ative compressibilities < can indeed arise, because 
the surface separating the two phases also contributes 
to the total energy. (From the occurrence of PS, it can 
be inferred that the phase boundary is not energetically 
favorable because the system would otherwise tend to 
maximize instead of minimizing it and form many small 
droplets or a mixed phase.) If this boundary has high 
enough energy and grows with doping, its contribution 
to the total energy leads to a negative compressibility. 

In the thermodynamic limit L — )■ oo, such surface 
terms become negligible compared to the bulk contri- 
butions, implying — > 0. Moreover, the boundary 
surface does not have to grow with doping at all: In the 
present one dimensional case, the phase boundary always 
consists of just two bonds connecting the two phases re- 
gardless of their size. For this reason, one could obtain 
« here even for finite systems. In addition, one 
has to keep in mind that numerical differentiation is no- 
toriously sensitive to even small numerical errors of the 
ground state energies entering Eq. 10}. 

A discontinuity of the density n as function of the 
chemical potential fi is equivalent to the limit of infi- 
nite compressibility k^^ — > 0, and is likewise taken as 
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FIG. 10: Grand canonical expectation value for the electron 
density (n) vs. chemical potential fi for J' — 0.01, L = 24 and 
(a) U — 0, (b) U = 10. The dashed lines indicate the limits of 
the (a) polaronic or (b) phase separated and the FM regimes, 
respectively. 



an indication for phase separation. We will here show 
that it can arise from independent polarons as well. 
In order to obtain from the DMRG calculations 

at fixed particle numbers N^, we set = N^/L, 

where minimizes the grand canonical expectation 
value (H — ^N) — E{Ne, L) — fiNg. The results for chains 
with L = 24, J' = 0.01, [/ = and [/ = 10, respectively, 
are shown in Fig. 1101 

We find a jump near n « 1 both for U ~ 10, were PS 
is observed, as well as for U ~ 0, where we see polarons. 
In the latter case, it can be accounted for by the inde- 
pendence of the polarons at low doping, i.e., adding each 
polaron costs the same energy, which is balanced by the 
chemical potential fj,* = epoi.-^ This also corresponds to 
our observation that the ground state energy per site as 
a function of filling (not shown) lies practically on an 
straight line near n — 1 with dE/dn — epoi. 

In conclusion, we have argued that a negative com- 
pressibility for a finite system is neither a necessary nor 
a sufficient condition for phase separation. It is also im- 
portant to note that a discontinuity of the density n as 
function of the chemical potential /i is a necessary condi- 
tion for phase separation, but is not a sufficient condition 
because other mechanisms such as independent polarons 
can induce a 'jump' of n{^) as well. 



IV. CONCLUSIONS 

We have presented an extensive numerical study of the 
one-dimensional ferromagnetic Kondo lattice model with 
quantum mechanical S — 3/2 core spins using the density 
matrix renormalization group method, treating a non- 
degenerate conduction band with and without on-site 
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Coulomb repulsion at low to intermediate doping. In par- 
ticular, we have explored the similarities with the anal- 
ogous model with classical core spins, where ferromag- 
netic polarons have been found to dominate over phase 
separation^ for parameters relevant to manganites. We 
have investigated whether the inclusion of quantum fluc- 
tuations leads to an attractive interaction between the 
polarons and thus to phase separation. We find that this 
is not the case: the polarons are in fact repulsive at short 
distances. In general, the results of the quantum model 
agree very well with the classical model and show that 
classical spins are indeed a very good approximation for 
manganites. 

Furthermore, we have investigated the influence of a 
relatively large local Coulomb repulsion U — 10, which 
reduces the effective AFM coupling J^s relevant when 
n < 1. For small t2g superexchange J' < 0.01, phase 
separation dominates over separated polarons if on-site 
Coulomb repulsion is taken into account, while polarons 



are favored upon increasing the superexchange parame- 
ter to J' = 0.02. We have summarized these findings in 
phase diagrams in the plane of doping and t2g superex- 
change J' for [/ = and U = 10. 

We have also argued that the observation of a (small) 
negative compressibility and of a jump in the density n as 
a function of the chemical potential /i are not compelling 
indicators of phase separation. While a discontinuity of 
n(/i) is a necessary condition for phase separation, it is 
not a sufhcient one: different mechanisms such as the for- 
mation of independent polarons can leave similar traces. 
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